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The aim of this short lecture course is to develop Selberg's trace formula for 
a compact hyperbolic surface M, and discuss some of its applications. The 
main motivation for our studies is quantum chaos: the Laplace-Beltrami oper- 
ator —A on the surface M represents the quantum Hamiltonian of a particle, 
whose classical dynamics is governed by the (strongly chaotic) geodesic flow 
on the unit tangent bundle of A4 . The trace formula is currently the only avail- 
able tool to analyze the fine structure of the spectrum of —A; no individual 
formulas for its eigenvalues are known. In the case of more general quan- 
tum systems, the role of Selberg's formula is taken over by the semiclassical 
Gutzwiller trace formula jllj , [7] . 

We begin by reviewing the trace formulas for the simplest compact man- 
ifolds, the circle §^ (Section [T]) and the sphere §^ (Section [21) . In both cases, 
the corresponding geodesic flow is integrable, and the trace formula is a con- 
sequence of the Poisson summation formula. In the remaining sections we 
shall discuss the following topics: the Laplacian on the hyperbolic plane and 
isometrics (Section [3]); Green's functions (Section [4|); Selberg's point pair in- 
variants (Section[5]); The ghost of the sphere (Section[6]); Linear operators on 
hyperbolic surfaces (Section [7]); A trace formula for hyperbolic cylinders and 
poles of the scattering matrix (Section [8]); Back to general hyperbolic surfaces 
(Section [9]); The spectrum of a compact surface, Selberg's pre-trace and trace 
formulas (Section I10|) : Heat kernel and Weyl's law (Section 111^: Density of 
closed geodesies f Section I12|): Trace of the resolvent (Section I13|) : Selberg's 
zeta function fSection ll4p : Suggestions for exercises and further reading (Sec- 
tion [TH]). 

Our main references are Hejhal's classic lecture notes [131 Chapters ONE 
and two], Balazs and Voros' excellent introduction [l, , and Cartier and Voros' 
nouvelle interpretation 6:. Section [T5l comprises a list of references for further 
reading. 

These notes are based on lectures given at the International School Quan- 
tum Chaos on Hyperbolic Manifolds, Schloss Reisensburg (Giinzburg, Ger- 
many), 4-11 October 2003. 
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1 Poisson summation 



The Poisson summation formula reads 



for any sufficiently nice test function : K ^ C. One may for instance take 
h e C^(M) with \h{p)\ < (1 + \p\y^^^ for some 6 > Q. (The notation x < y 
means here there exists a constant C > such that x < Cy.) Then both sums 
in ([T]) converge absolutely. ([T]) is proved by expanding the periodic function 

f{p)=Y.h{p + m) (2) 

in its Fourier series, and then setting p = 0. 

The Poisson summation formula is our first example of a trace formula: 
The eigenvalues of the positive Laplacian —A — —-^ on the circle § of 
length 27r are where m = 0, ±1, ±2, . . ., with corresponding eigenfunctions 
fm{x) = (27r)~^/^e'™^. Consider the linear operator L acting on 27r-periodic 
functions by 

[Lf]{x) / \ix,y)fiy)dy (3) 
Jo 

with kernel 

Hx,y) = E K'm)fm{x)^^{y). (4) 



m£2 



Then 



Liprn = h{m)ip„i (5) 
and hence the Poisson summation formula says that 

The right hand side in turn has a geometric interpretation as a sum over the 
periodic orbits of the geodesic flow on whose lengths are 27r|n|, n E Z. 

An important example for a linear operator of the above type is the re- 
solvent of the Laplacian, {A + with Imp < 0. The corresponding test 
function is h{p') — {p^ — p'^)^^ ■ Poisson summation yields in this case 

r „ — 27ri|ji|p' 



and by shifting the contour to — loo and collecting the residue at p' = p, 
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Imp' 



Rep' 



Q 



we find 



2 m2)-i 



-27ri|n|/9 



(8) 

The right hand side resembles the geometric series expansion of cotz for 
Im < 0, 



2ie cos 2; n,,. _ 

cot z = — = i(l + e ^'"^) e 

n=0 



2inz 



1 - 



ne2 



Hence 



which can also be written in the form 



2\-l 



mez 



1 1 

+ 



p — m p + m 



■ cot(7rp), 



= TT cot(7rp), 



(9) 
(10) 

(11) 



The above h is an example of a test function with particularly useful analytical 
properties. More generally, suppose 

(i) h is analytic for | Im p| < cr for some cr > 0; 

(ii) \h{p)\ <C (1 + I Rep|)~^~'' for some 6 > 0, uniformly for all p in the strip 
I ImpI < a. 



Theorem 1. If h satisfies (i), (ii), then 

Yl '^(™) ^^P'> cot(7rp) dp 

mez ' "'<^= 

where the path of integration C= is 



(12) 





Imp 




Rep 


4 






► 
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Proof. The Poisson summation formula ([T]) may be written in the form 

J2 h{m) ^\Y.I [Mp) + H-P)] e-'^'l-l" dp. (13) 



We shift the contour of the integral to J^^^^^. The geometric series expansion 
of cot z in ([9]) converges absolutely, uniformly for all z with fixed negative 
imaginary part. We may therefore exchange summation and integration, 

2^ /'OC — 1(7 
^^^77 -/-OO-lfT 



We conclude with the observation that 



— 00-|-l(T 



— i h{~p) cot(7rp) dp = - I h{p) cot(7rp) dp (15) 

since cot z is odd. □ 

Remark 1. This theorem can of course also be proved by shifting the lower 
contour in (fT2l) across the poles to the upper contour, and evaluating the 
corresponding residues. 



2 A trace formula for the sphere 

The Laplacian on the sphere is given by 

sm oO \ oO J sm 9 0(p 

where G [0,7r), e [0,27r) are the standard spherical coordinates. The 
eigenvalue problem 

{A + A)/ = (17) 

is solved by the spherical harmonics / = YJ'" for integers I — 0,1,2,..., 
m = 0, ±1, ±2, . . . , ±1, where 

1/2 

P;"(cos6i)e'™'^ (18) 

and P™ denotes the associated Legendre function of the first kind. The eigen- 
value corresponding to is A = + and hence appears with multiplicity 
2Z + 1 . Let us label all eigenvalues (counting multiplicity) by 

Aq < Ai < A2 <...-> 00, (19) 



(2/ + 1) {I -my. 
4tt {l + m)\ 
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and set pj = y Aj + i > 0. For any even test function h e C^(R) with the 

bound \h[p)\ ^ (1 + | Rcp|)^^^* for some (5 > (assume this bound also holds 
for the first and second derivative) we have 

oo oc 

Y.Hp,)=Y,{2l + l)h{l + \) (20) 

3=0 1=0 



E (21) 

1 — — 00 

oo „ 

E / |? + i|M^ + 5)e'"''"rf/ (22) 

/>oo 

2E(-1)" / ph{p) cos{2Trnp)dp, 



(23) 



in view of the Poisson summation formula. We used the test function |p|ft,(p) 
which is not continuously differentiable at p = 0. This is not a problem, since 
we check that (using integration by parts twice) 

p h{p) cos(27rnp) dp = 0{n-'^) (24) 

hence all sums are absolutely convergent. With Area(§^) = 47r, it is suggestive 
to write the trace formula ([23]) for the sphere in the form 



EMft) = ^^^^ / |p|Mp)rfp + E(-l)" / |p|Mp)e''^'"^dp. (25) 

As in the trace formula for the circle, the sum on the right hand side may 
again be viewed as a sum over the closed geodesies of the sphere which, of 
course, all have lengths 27r|7i|. The factor (—1)" accounts for the number of 
conjugate points traversed by the corresponding orbit. 

The sum in (j23p resembles the geometric series expansion for tan z for 
Im z < 0, 

tanz = -cot(z + 7r/2) ^ -i E(-l)"e~^'l"l^ (26) 

As remarked earlier, the sum converges uniformly for all z with fixed Im z < 0. 
We have in fact the uniform bound 

E (-l)"e"^'l"l^ < l + 2Ee^"^"'^ < 1 + 2 / e^^^^^da;^: (27) 

nGZ n=l •'^ ™^ 

which holds for all z with Imz < 0. 

Let us use this identity to rewrite the trace formula. Assume h satisfies 
the following hypotheses. 
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(i) h is analytic for | Im p\ < c for some u > 0; 

(ii) h is even, i.e., h{~p) — h{p); 

(iii) \h{p)\ <C (1 + I Re p\)^'^~^ for some 6 > 0, uniformly for all p in the strip 
I Im p\ < cr. 



Theorem 2. If h satisfies (i), (ii), (iii), then 

^ n C V 



(28) 



wiii/i i/ie pat/i of integration 





Imp 




a 

















Proof. We express (1^51) in the form 

poo poo 

^(-1)"/ pMp)e^-'l"l''dp + ^(-l)" / p%)e-2-l"l''dp. (29) 
nez •^'^ nez •^'^ 

which equals 

/O /'OO 
p/.(p)e-2H"|Prf^ + ^(_l)« / p/i(p)e'2H«|Prf^. (30) 
-00 Jo 



Let us first consider the second integral. We change the path of integration to 
Ci: 



^--^l 




Imp 




a 






Rep 




!/ 






C, 



Due to the uniform bound ((27|) and 

/ \ph{p)[l - (27rImp)-Vp| < 



(31) 



we may exchange integration and summation, and hence the second integral 
in ([50)1 evaluates to 
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^(-1)" / p/i(p)e~2^'l"l''dp = i / hip)ptaTi{7Tp)dp (32) 



The first integral in (jSO]) is analogous, we have 



-^(-1)"/" p/i(p)e-2'^'l"l''dp== -1 / /i(p)ptan(7rp)dp (33) 

= 1 / h{p) ptaii{TTp) dp. (34) 

The final result is obtained by reflecting these paths at the origin, using the 
fact that h is even. □ 

Remark 2. The poles of tan z and corresponding residues can be easily worked 
out from pT|) . 



7rtan(7rp) = — Trcot [7r(p+ ^ 

OC r 



i=-c 



^— — oo 



1 



p-il-^) P+il + k) 

1 1 
P+il+'^)^ p-{l + k) 



(35) 
(36) 

(37) 



(the sum has not been reordered, we have simply shifted the bracket) 

^ " 1 1 " 



E 

/=0 



p+(? + i) p-(/ + i) 



(38) 



Note that the extra factor p in the integral (|28|) . as compared to Theorem [TJ 
yields the multiplicity of the eigenvalues of the sphere. 



3 The hyperbolic plane 

In this section we briefly summarize some basic features of hyperbolic geom- 
etry; for a detailed discussion see Buser's lecture notes [5]. 

The hyperbolic plane may be abstractly defined as the simply con- 
nected two-dimensional Riemannian manifold with Gaussian curvature —1. 
Let us introduce three convenient coordinate systems for H'^: the Poincare 
disk D = {z : \z\ < 1}, the upper half plane Sj = {z : Imz > 0} and polar co- 
ordinates (r, (f)) e M>o X [0, 27r). In these parametrizations, the line element ds, 
volume element dfi and the Riemannian distance d{z, z') between two points 
z, z' G are as follows. 
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(1 — x'^ — 



coshd(2;, z') 



4 dx dy 



2\z-z'\ 



dx'' + dy'' 

jj2 



(1 _ a;2 _ ^2)2 



(1 - \z\ni - \z'\) 



Tvr 



dxdy 



1 



21mzlmz' 



polar dr + sinh r 



sinh T dr ( 



coshr [for z — (r, 



z 



(0,0)] 



The group of isometries of H^, denoted by Isom(EI^), is the group of smooth 
coordinate transformations which leave the Riemannian metric invariant. The 
group of orientation preserving isometries is called Isom"'"(H^). We define the 
length of an isometry g £ Isom(EI^) by 



= ^{9) = ^fi^ d{gz, z). 



(39) 



Those g G Isom^(]HI^) for which £ > are called hyperbolic. In the half plane 
model, Isom'''(lI^) acts by fractional linear transformations, 



g:Sj^ 9), 



gz := 



az + b 
cz + d^ 



c d 



G SL(2, 



(40) 



(we only consider orientation-preserving isometries here). We may therefore 

identify g with a matrix in SL(2.]R), where the matrices g and —g obviously 
correspond to the same fractional linear transformation. Isom^(IHI^) may thus 
be identified with the group PSL(2,]R) = SL(2,M)/{±1}. In this representa- 
tion, 

2cosh(^g/2) = max{|tr.g|,2}, (41) 
since every matrix g G SL(2,M) is conjugate to one of the following three. 



1 b 
1 



a 
a-i 



cos(6l/2) sin(6'/2)\ 
-sin(6l/2) cos(6'/2)j ' 



(42) 



with 6 G R, a G R>o, G [0, 2tt). 

The Laplace- Beltrami operator (or Laplacian for short) zi of a smooth 
Riemannian manifold with metric 



ds^ = ^y^gjkdx^dx^ 



(43) 



is given by the formula 



(44) 



where g^'^ are the matrix coefficients of the the inverse of the matrix {gjk)^ 
and g = \Act{gjk)\- In the above coordinate systems for the Laplacian 
takes the following form. 
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A 


D 


{l-x' - 
4 






dy^ ) 








( d'^ 


f 






polar 


1 d 






1 




sinh r dr 


sinh^ T 


502 



One of the important properties of the Laplacian is that it commutes with 
every isometry g g Isom(H2). That is, 

ATg ^ TgA V.g e Isom(H2). (45) 

where the left translation operator Tg acting on functions / on H2 is defined 

by 

[Tgf]{z) = f{g-'z) (46) 

with g G Isom(H2). Even though (I45p is an intrinsic property of the Lapla- 
cian and is directly related to the invariance of the Riemannian metric under 
isometrics, it is a useful exercise to verify (|^5]) explicitly. To this end note that 
every isometry may be represented as a product of fractional linear transfor- 
mations of the form z i— > az (a > 0), z i-^ z -I- 6 (6 € M), z i— > — 1/z, z i— > — z. 
It is therefore sufficient to check p5)) only for these four transformations. 



4 Green's functions 

The Green's function G(z,w;A) corresponding to the differential equation 
[A + A)/(z) = is formally defined as the integral kernel of the resolvent 
{A + A)^^, i.e., by the equation 

{A + \)-^f(z) ^ J G{z, w- A) f{w) dfiiw) (47) 

for a suitable class of test functions /. A more precise characterization is as 
follows: 

(Gl) G( • , w; A) e C°°(H2 - {w}) for every fixed w; 
(G2) {A + X}G{z, w; A) = S{z, w) for every fixed w; 

(G3) as a function of (z, w), G{z, w; A) depends on the distance d(z, w) only; 
(G4) G(z, w; A) ^ as d(z, w) oo. 

Here S{z, w) is the Dirac distribution at w with respect to the measure dfi. It 
is defined by the properties that 

(Dl) S{z,w)dii{z) is a probability measure on IHl2; 
(D2) /jj, /(z) S{z, w) dii{z) = f{w) for all / € C(M2). 
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E.g., in the disk coordinates z — x + ly, w — u + iv £ Ti we then have 

(1 _ a;2 _ y2\2 

S{z, iv) = S{x — u) 5{y — v) 



(48) 



where S{x) is the usual Dirac distribution with respect to Lebesgue measure 
on M. In polar coordinates, where w is taken as the origin, t = d{z, w), and 



S{z, w) 



Sir) 



2tt sinh r 

Property (G2) therefore says that (A + X)G{z, w; A) = for z w, and 



{A + X)G{z,w;X)dn{z) ^1 Ve > 0. 



(49) 



(50) 



d{z.w)<e 



In view of (G3), there is a function / G C°°(M>o) such that /(r) — G{z, w; A). 
Then 



1 = 



d{z,'w)<e 



{A + X)G{z,w;X)d^{z) 



ZTT — smhr — + Asnih 
Jo \dT V drj 

27rsinhe/'(e) + 27rA / sinh T/(r) dr. 

"'0 



r /(r) dT 



(51) 
(52) 
(53) 



Taylor expansion around e = yields /'(e) — l/(27re) + 0(1) and thus /(e) 
(l/27r) loge + 0(1) as e ^ 0. Equation (G2) is therefore equivalent to 



{A + X)G{z,w;X) = 0, 



d{z, w) > 0, 



[G(z,w;A) = (l/27r)logd(z,w) +0(1), d{z,w) 0. 

Proposition 1. If p E C with Imp < 1/2, and X = p'^ + j, then 

1 



(54) 



G{z,w;X)=-—Q_ 



coshd(z, w)) (55) 

satisfies (G1)-(G4), where Qu is the Legendre function of the second kind. 
Proof With /(r) = G(z, w; A), becomes 



1 



sinh T dr \ 



( sinhr— ) + A 



fir) = 0. 



(56) 



Set r = coshr, /(coshr) — fir), and A — —v{v + 1), to obtain Legendre's 
differential equation 



^ ^ dr^ 



9, d 

2r—+viv + l) 
dr 



fir) = 0, 



(57) 
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whose solutions are the associated Legendre functions Pv{r) and Qi,{r). 
has the integral representation 



1 r 

Q_.+.,(coshr)^-^| 



V cosh t — cosh T 



dt (58) 



which converges absolutely for Imp < 1/2 and r > 0. From this it is evident 
that Q_ i_^j^(coshr) ^ as r ^ oo (see also Lemma [1] below), so (G4) holds. 
For i — > (p fixed) it has the asymptotics required in (|54p . cf. the well known 
relation 

Q_i+,p(coshT) ~ -(log(r/2)+7 + ^.(i +ip)) (59) 

where 7 is Euler's constant and the logarithmic derivative of Euler's F 
function. □ 

Lemma 1. Given any e > there is a constant Cg > such that 

|Q_i+,,(coshT)| < C,(l + I logr|)e-(i'"''-5+^) (60) 

uniformly for all t > and p G C with Imp < i — e. 
Proof. From the integral representation (|58p we infer 

gt Imp 



1 

|Q_i+.p(coshr)|<^| 



V cosh t — cosh r 



(61) 



smce 



1 1 f 0^2' 

< 1 gr(Imp-2+^) / =dt (62) 

V2 Jr vcosht — coshr 



g*(Imp-|+«) < gT(Imp-i+e) ^gg-j 



for all t > T, if Im p— i + e<Oas assumed. The remaining integral 



e*(|-e) 

= (64) 
T V cosh t — cosh T 

has a logr singularity at t = and is otherwise uniformly bounded for all 

T — > 00. □ 

Remark 3. This is only a crude upper bound, but sufficient for our purposes. 
To highlight the p dependence, we shall use 

Gp{z,w) = -^g_i+jp(coshd(z, w)) (65) 



instead of G(z, w; A). 
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Lemma 2. Suppose / : — > C with \f{z)\ < Ae"'^''^'°\ with constants 
A,a>0. Then the integral 

I Gp{z,w)f{w)df,iw) (66) 

converges absolutely, uniformly in Rep, provided Imp < — (a + i). The con- 
vergence is also uniform in z in compact sets in . 

Proof. Note that |/(w)| < yle"''(°^"') < ^e"''(°'^)e"''(^''"). In polar coordinates 
(take w as the origin) r = d{z,w), and dfi = smhr dr d(j). In view of Lemma 
m the integral (|66)) is bounded by 

poo ^ 

<e / llogrl e-(2-i'"''-^-")-sinhrdT (67) 
Jo 

which converges under the hypothesis on Im p. □ 



5 Selberg's point-pair invariants 

Let H he a, subgroup of Isom(]HI^). An i/-point-pair invariant fc : x ^ C 
is defined by the relations 

(Kl) k{gz,gw) ~ k{z,w) for all g E H, z,w H^; 
(K2) k{w, z) = fc(z, w) for aU z, w e H^. 

Here we will only consider point-pair invariants which are functions of the 
distance between z,w, such as the Green's function Gp{z,w) studied in the 
previous section. Hence H = Isom(H^) in this case. We sometimes use the 
notation fc(T) = k{z,w) with t = d{z,w). Let us consider 

1 

k{z,w) ~ — / Gp{z,w)ph{p) dp (68) 

""l J-oo 

where the test function h : C ^ C satisfies the following conditions. 

(HI) h is analytic for | Imp| < cr for some a > 1/2; 
(H2) h is even, i.e., h{—p) = h{p); 

(H3) \h{p)\ < (1 + iRepD^^^'^ for some fixed (5 > 0, uniformly for aU p in 
the strip jlmpj < a. 

For technical reasons we will sometimes use the stronger hypothesis 

(H3*) \h{p)\ <jv (1 + I Rcp|)"^ for any fixed iV > 1, uniformly for aU p in 
the strip | Impj < a. 
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The Fourier transform of h is 

9{t)^^ f h{p)e-^'"dp. (69) 

JR 

With the integral representation ([58|) one immediately finds 
1 f°° a'(t) 

k{z,w)^ y= =dt, T = d{z,w). (70) 

Try 2 Jt V cosh t — cosh r 

The analyticity of h and (H3) imply that g and its first derivative {all deriva- 
tives provided (H3*) holds) are exponentially decaying for \t\ oo. To see 
this, consider 

5^^^W = 7^ /(-vrMp)e-'^'rfp (71) 



27r 
1 

2^ 
1 

2^ 



{-iprhip)e-^'''dp (72) 
-i{p-i(j)Yh(p-ia)e-'P*dp. (73) 



Since, due to (H3*) 



wc find 



/ \ip - lay h{p-ia)\ dp <oo (74) 



|ff('')(i)|«.e-'^l*l. (75) 
The point-pair invariant k(z,w) gives rise to the linear operator L defined 



by 



[Lf]{z) / k{z,w)fiw)df,{w). (76) 

JH2 



Proposition 2. Suppose f G C^(]HI^) is a solution of (Z\ + + i)/ = with 
I Imp| < cr anrf |/(z)| < Ae°"^^''^°\ with constants A> 0,0 < a < a-^. Then, 
for h satisfying (HI), (H2), (H3), 

Lf = Kp)f. (77) 

Proof. We have 

(78) 

Gp\z,w)p'h{p')dp' \f{w)diji{w) (79) 



/ 


/fc(; 
















/ 




'H2 



— oo — icr 



Gp. (z, w;)p' /i(p') dp' ^/(u.)dAi(w;) (80) 
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where we have shifted the contour of integration by <j. Then 

/•oo — icr 



1 

TTl 



— oo — ICT 



Gpi (z, w)f{w)dfi{'w) 'fp' h{p') dp (81) 



since the inner integral converges absolutely, uniformly in Re p, see Lemma [2j 
We have 

GAz,yj)fMdp{w) = {A + p'' + l)-V(^) = ip'' - P^)-'f{z) (82) 
and thus 



\Lm = -^!{z) I 4?^*'. (83) 



I"! J-oo-ia P - P^ 

This integral converges absolutely, cf. (H3), and is easily calculated. We shift 
the contour from ~icr to +icr and collect residues, so that 



1 f°°-'^ 2p'h{p') , , , , 1 /■°°+'" 2p'h{p') ^ , 

^''1 J-00-l(T P — P ^''1 J-OO + IO- P ~ P 

Since h is even, the integral on the right hand side equals the negative of the 
left hand side, and thus 

1 p-'^ 2p'h{p') 

which concludes the proof. □ 

It is useful to define the auxiliary functions ^, Q : M>o ^ C by the relations 
<?>(2(coshr-- 1)) = fc(r) = fc(z,w), t = d{z,w), (86) 

and 

g(2(coshi- 1)) =5(t). (87) 
Lemma 3. The following statements are equivalent. 

(i) h satisfies (HI), (H2), (H3*). 

(ii) g e C°°(R>o) with 

|Q('^)(77)| (1 + 7])"""" V?7>0. (88) 

Proof. Clearly g € C°°(R) if and only if Q e C°°(R>o) (this is obvious for 
t ^ 0; the problem at i = can be resolved by expanding in Taylor series). 
In view of ([75)1 . the bound ([55)1 is evident for v — Q. The i^th derivative of g 
is of the form 

g(-)(i) =^aj,e^'l*lQ(j')(2(cosht-l)) (l + 0(e^l*l)) (89) 
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with suitable coefficients a^v Hence, by induction on v, 

e^tl |QM(2(cosh<- 1))| <^ \9^''\t)\ + | ^a^^e^'l*' Q^^'^ (2(cosht - 1)) | 

(90) 

W| +^eJl*le(-"-^)l*l (91) 

«. e-"l*l. (92) 

This proves (i) (ii) . Conversely, (|88p implies via ([M]) the exponential decay 
of 5, which proves (HI). (H3*) follows from g e C°°(M). □ 

The integral transform ([70)1 reads in terms of the functions Q,??, 

In order to specify Q uniquely for a given <P, we always assume in the following 
that (3(7?) — » for 77 oo. 

Lemma 4. Consider the following conditions. 

(i) g e C°° (R>o) and | gf") (r/) | (l + ry) """" . 

(ii) e C°°(M>o) and |<Z'('^)(0| «. (1 + e)"""""'^'^'. 

T/ien (i) implies (ii) /or on?/ /iierf e > 0, and (ii) implies (i) /or any /ixed 
e < 0. 

Proof. The I'th derivative of (p is 

di" Jo 

dr^. (96) 

Therefore (i) imphes G C°°(R>o). Furthermore, from ([^5)) . 



Vv 

1 p g(''+i)(r?) 



Jo 

Jo Vl 
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for e > small enough. The last integral converges for any e > 0. 

The implication (ii) (i) follows analogously from the inversion formula 

Qi^)^r^=dc (99) 

To show that (l99l) is indeed consistent with (l93t. write (l99l) in the form 



and thus 



/oo 
<P{v + e)d^ (100) 
-oo 



/OO 
'P'{v + e)dt (101) 
-OO 

The right hand side of ^ is 



_1 rO^^^^^^i f Q'(^ + ^2)^^ (102) 

= --( f <P'{( + T^l + 7^l)dr^idm, (103) 
Jr Jr 

where we have used (jlOip in the last step. This equals of course 

) POO 

■P' {£, + r^)rdr = - <P' + r)dr = <P{^) (104) 



which yields the left hand side of (|93p . □ 

Proposition 3. The following statements are equivalent. 

(i) h satisfies (HI), (H2), (H3*). 

(ii) k{z,w) is in C°°(]HI^ x H^) with the bound on the vth derivative, 

Ifc^'H^)! e-("+i/2-^)" Vt>0, (105) 
for any fixed e > . 

Proof. In view of Lemmas [3] and SI the statement (i) is equivalent to the condi- 
tion for statement (ii) in LemmaS) Since fc(z, w) = fc(r) = ^(2(coshT— 1)) , 
the proof is exactly the same as that of Lemma |3] with g(t) replaced by fc(r), 
and Q(2(cosht - 1)) by ^(2(coshT - 1)). □ 



6 The ghost of the sphere 



Note that for z = w, the kernel k{z, w) has a finite value, unlike the logarithmic 
divergence of the Green's function Gp{z, w). In fact, 
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fc(z, z) = — - 



g'it) 



■kV2 Jo \J cosh t—\ 



dt 



1 



9'{t) 



27r Jo sinh(i/2) 



dt 



1 
1 

4^ 



sin(pt) 
sinh(t/2) 
sin(pt) 



ft-(p) pdp j-dt 
dt >h{p) pdp 



/o sinh(t/2) 

where changing the order of integration is justified, since 







sin(pt) 



sinh(t/2) 



dt < \p\ 



sinh(i/2) 



dt < \p\ 



17 
(106) 
(107) 
(108) 
(109) 

(110) 



and \h{p) \ ^ (1 + |p|) assuming (H3*). We use the geometric series expan- 
sion 

1 2e-*/2 



sinh(</2) 1 - e-* 
so for I ImpI < 1/2 

sm{pt) 



2^exp [- + t], 



1=0 



sinh(i/2) 



dt^^Yl 



1=0 



ip-(z + i) ip + (; + i) 



-7ritan(7ri/3) = 7rtanh(7rp), 



compare 



We conclude 



1 f°° 

kiz,z) — — / h{p) ta,nh{TT p) pdp. 
47r J_^ 



(111) 

(112) 
(113) 

(114) 



Let us conclude this section by noting that the logarithmic divergence of 
the Green's function is independent of p, see ([54]) . It may therefore be removed 
by using instead 

Gp{z,w)~GpAz,w) (115) 
where p, 7^ p is a fixed constant in C with | Im p^\ < 1/2. We then have from 



lim [Gp{z, w) — Gp^ (z, w)] — 



1 



2'K\/2 Jt V cosh t — cosh r 



dt 



1 



00 _ipt 



Att 7o sinh(i/2) 

-E 



27ri 



p-i(; + i) p*-i(Z + i) 



(116) 
(117) 

7 

(118) 
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where we have used the geometric series expansion (|llip as above. The last 
sum clearly converges since 

^-TT 77-TT-0(r2)- (119) 

Remark 4- In analogy with the trace formula for the sphere, we may view the 
geometric series expansion, 

tanh(7rp) = ^(-l)"e-2^l"l^ (120) 

cf. ((26)) . as a sum over closed orbits on the sphere, but now with imaginary 
action. These orbits have an interpretation as tunneling (or ghost) orbits. 



7 Hyperbolic surfaces 

Let Al be a smooth Riemannian surface (finite or infinite) of constant negative 
curvature which can be represented as the quotient r'\H^, where r' is a strictly 
hyperbolic Fuchsian group (i.e., all elements 7 G -T — {1} have > 0). The 
space of square integrable functions on Ai may therefore be identified with 
the space of measurable functions / : — )■ C satisfying the properties 

T^f = f V7er (121) 

and 

||/||2:= / \f\'d^i<^ (122) 

where is the translation operator defined in (|46p and J-'p is any fundamental 
domain of in H^. We denote this space by L^(r\H2). The inner product 

(/i,/2)= / /JadM (123) 

makes L^{r\lP) a Hilbert space. Similarly, we may identify C°°(I^\IHI^) with 
the space of functions / G C°°{M'^) satisfying (|12ip (note that more care has 
to be taken here when F contains elliptic elements). Since A commutes with 

T^, it maps c°°ir\M^) c°°{r\m^). 

To study the spectrum of the Laplacian on /^\H^, let us consider the linear 
operator L of functions on P\H^, 

[Lf]{z):^ f kr{z,w)f{w)diJL{w) (124) 
Jr\vp 

with kernel 

= ^fc(7z,7«) (125) 

with the point-pair invariant k as defined in (|68p. The convergence of the sum 
is guaranteed by the following lemma, cf. Proposition 0] below. 
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Lemma 5. For every 5 > Q, there is a C's > such that 



J2 e-(i+^)'^('^^."') < Cs (126) 
■yer 



for all {z,w) e x 



Proof. Place a disk ^^{r) = {z' E H : d{'-^z,z') < r} around every point 
z-y — 7z, and denote the area of Ti^ir) by Area(r). Then 

< f''^'^ I e-(i+^)'^(^»dM(z')- (127) 
- Area(r) 

If 

r < - min £^ (128) 
2 7er-{i} ^ ^ ' 

the disks f-y(r) do not overlap. (Note that min^£/-_{ij > since F is 
strictly hyperbolic and acts properly discontinuously on H^.) Therefore 



y g-(i+5)d(7..») < /■ e^(i+*')''(-''-)dM(z') (129) 

^ - Area(r) 

This integral converges for any (5 > 0. □ 

Proposition 4. // /i satisfies (HI), (H2), (H3*), t/ieri f/ie kernel kr{z,w) is 
in C°"{r\m^ X AH2), with kr{z,w) = kr{w,z). 

Proof. PropositionOand LemmaOshow that the sum over k{^z, w) converges 
absolutely and uniformly (take 5 = <t — 1/2 — e > 0). The same holds for sums 
over any derivative of k{'jz,'w). Hence kr{z,w) is in C°°(IEII^ x H^). To prove 
invariance under note that 

krijz^w) = k{j'jz,w) — k{'-f'z,w) = kr{z,w). (131) 
7'e-r 7'e-r 

Thus kr{z,w) is a function on F\H^ with respect to the first argument. 
Secondly 

kr{z,w) = k{'j'z,w) — k{w,j'z) 
I'er 7'er 

= Y Hl'~^w,j'^^-/'z) ^ Y k{-/'^^w,z) ^ kr{w,z), (132) 
7'e-r 7'6_r 

which proves symmetry. Both relations imply immediately kr{z,"fw) — 
kr(z,w). □ 
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Proposition 5. Suppose f e C'^{r\E.'^) is a solution of {A + + i)/ = 
with |Imp| < cr and \ f{z)\ < Ae°"^^''^°\ with constants A> 0,0 < a < a ~ 
Then, for h satisfying (HI), (H2), (H3), 

Lf = h{p)f. (133) 

Proof. Note that 

/ kriz,w)f{w)dij{w) — / k{z,w)f{w)dp{w) (134) 
and recall Proposition [21 □ 



8 A trace formula for hyperbolic cylinders 

The simplest non-trivial example of a hyperbolic surface is a hyperbolic cylin- 
der. To construct one, fix some 7 S Isom^(]HI^) of length ^ — £{-f) > 0, and 
set r — ^, where 3 is the discrete subgroup generated by 7, i.e., 

3 = {7" : neZ}. (135) 

We may represent 3\IHI^ in halfplane coordinates, which are chosen in such a 
way that 

7 = (^0' . (136) 

A fundamental domain for the action of 7 on io, z 1— > e^z, is given by 

{z e : 1 < y < e^}. (137) 

It is therefore evident that has infinite volume. A more convenient set 

of parameters for the cylinder are the coordinates (s,u) G M^, with 

X = ue", y = e^ (138) 

where the volume element reads now 

dp — ds du. (139) 

In these coordinates, the action of 7 is (s, m) i— > (s + u), and hence a funda- 
mental domain is 

= {(s, u) e : < s < ^}. (140) 

Note that 

coshd(7"0,z) = 1+ I— -^_L = l + 2sinh2(n£/2)(l + ii2) (141) 



Selberg's trace formula: an introduction 21 

and hence 

oc 

fc3(z,z) = fc(z,2) + ^fc(7"z,z) (142) 

oo 

= fc(z,z) + 2^/c(7"z,z) (143) 

n=l 
oo 

= fc(z,z) + 2^<?(4sinh^(n^/2)(l + M2)). (144) 

n=l 

From this we can easily work out a trace formula for the hyperbolic cylinder: 
Proposition 6. If h satisfies (HI), (H2), (H3), then 

( [fc3(z,z)-M^,z)]d/i = f;-M^. (145) 
J^\v? ^1 smh(n^/2) 



Proof. We have 



${Asmh^{ni/2){l + u'^))dsdu (146) 



~ 2sinh(n^/2) 

and with pOO)) . 

I 



j <?>(4sinh2(n^/2)+C^)fiC (147) 



Q(4sinh^(n^/2)) (148) 



2sinh(n£/2) 

which yields the right hand side of (fTl5l) . cf. (jST]). □ 
Proposition 7. ///i satisfies (HI), (H2), (H3), t/ien 

[/c3(z, z) — /c(z, z)](i/i = / h[p)n-^[p)dp (150) 
3\H2 7k 

w/iere 

I °° 

"3(P) = - E {°^P [("^ + I + 1/^) ^] - 1}"' (151) 

m=0 

is a meromorphic function in C wit/i simple poles at the points 

Pi^m = ^ + i{m+^) , j/eZ, m = 0,1,2,..., (152) 



and residues reSp^^ 713 = l/(7ri). 
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Proof. The geometric series expansion of 1/ sinh (jllip yields 

S Mw, ^ ; /. S S »p [- (■» + H w M.) * (153) 

and using again the geometric series, this time for the sum over n, 

oo 

^ exp [- (m + i + ip) n£] = {l - exp [- (m + i + ip) £] - 1 (154) 

n=l 

= {exp [(m+ i +ip)^] - (155) 
This proves the formula for 713 (p) . Near each pole p^m we have 

n^ip) ^ ^ {exp [(m + i + ip) ^] - l}"' (156) 



1 1 

TTl 



p- [(27r/£)i/ + i(m+i)] 
and so reSp^^ 71.3 = l/(7ri). □ 



(157) 



The poles of 713 (p) are called the scattering poles of the hyperbolic cylinder. 
A useful formula for 77,3 is 

"3(P) = :r E ■ ^""^P < V2, (158) 

ZTT smh(M/2) 

which follows immediately from the above proof. Furthermore, by shifting the 
path of integration to —100, we have the identity 

[ ^?^dp' = -n3ip), lmp<0. (159) 
Jr p2 _ p'^ p 



9 Back to general hyperbolic surfaces 

Let us now show that the kernel kr{z,w) of a general hyperbolic surface 
r\lP (with r strictly hyperbolic) can be written as a superposition of kernels 
corresponding to hyperbolic cylinders. 

Define the conjugacy class of any element 7 g _r as 

{7} := {7 e : 7 = g'^g^^ for some g e F}. (160) 

Clearly the length £^ is the same for all elements in one conjugacy class. The 
centralizer of 7 is 

3^:={ger: gj^ ^g}. (161) 



Selberg's trace formula: an introduction 23 

Lemma 6. If j F is hyperbolic, then the centralizer is the infinite cyclic 
subgroup 

3^ = {7" : ^ e Z}, (162) 
where "/^ € F is unique element such that 7™ = 7 for some m G N and there 
is no ^ F such that 7" = 7* for any n G N, n > 1. 

(The element 7* is called the primitive of 7.) 
Proof. 7 G PSL(2,R) ~ Isom^(H^) is conjugate to a diagonal matrix 





e-^/2 
with i-y > 0. The equation 

\ fab\ fa b\ fe^-'/'^ 



6-^/2/ \c d ~ \cd \ e-^-/2 



(163) 



(164) 



has the only solution b — c — Q, a — d^^ . Hence the centralizer is a diagonal 
subgroup of (a conjugate of) F . Since F is discrete, the centralizer must be 
discrete, which forces it to be cyclic. □ 

Remark 5. If 7 is hyperbolic and F strictly hyperbolic, then {7} ^ {7"} for 
all n 7^ 1. Furthermore, the centralizers of 7 and 7" coincide. 

The sum in (|125p can now be expressed as 

^k{-fz,w) = k{z,w) +^ ^ k{g'^-fgz,w) (165) 

= k{z,w)+^ ^ k{'ygz,gw) (166) 

where the respective first sums run over a set H of hyperbolic elements, which 
contains one representative for each conjugacy class {7}. We may replace this 
sum by a sum over primitive elements. If we denote by C H the subset of 
primitive elements, (jl66p equals 

= k{z,w) + ^ ^ k{jgz,gw) (167) 

OC 

= k{z,w)+Y, E E^(^V,5«^) (168) 

and hence, finally, 

kr{z,w) - k{z,w) = ^ E {hyigz, 9^) - k{gz, gw)} ; (169) 

recall that 

fc3^(z,u;) = E^(^"^'^)- (170) 

riGZ 
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10 The spectrum of a compact surface 

It is well known that for any compact Ricmannian manifold, —A has positive 
discrete spectrum, i.e., 

= Ao < Ai < A2 < .. . ^ oo, (171) 

with corresponding eigenfunctions (po = const, ipi, ip2, • ■ ■ G C°°(r'\]HI^), which 
satisfy 

(Z\ + A,)(^, =0 (172) 

and form an orthonormal basis of h^{r\lP). Furthermore, since A is real- 
symmetric, the ipj can be chosen to be real-valued. We furthermore define 

Pj = \A7~I. -7r/2<argpj <^/2. (173) 

If / e C^(r\H2), the expansion 

/(z)=^c,^,(z), c, = (/,^,), (174) 
j 

converges absolutely, uniformly for all 2; G H^. This follows from general spec- 
tral theoretic arguments, compare [131 p. 3 and Chapter three]. 

Proposition 8. If h satisfies (HI), (H2), (113), then 

Ltpj ^ h{pj)tpj . (175) 

Proof. Apply Proposition [5l Each eigenfunction tpj is bounded so a — 0. 
Furthermore, by the positivity of —A, we have | Imp| < 1/2 < a. □ 

Proposition 9. If h satisfies (HI), (H2), (H3*), then 

OQ 

kr{z,w) =:^^h{pj)ipj{z)lpj{w), (176) 
j=o 

which converges absolutely, uniformly in z,w E H^. 

Proof. The spectral expansion (|174|) of fcr(z, w) as a function of z yields 

00 

kr{z,w) ^^Cjipjiz), (177) 
i=o 

with 

Cj = / kr{z, w)Tpj{z) d^i{z) = [Lifij]{w) = h{pj)lpj{w) = h{pj)Tp^{w). 

(178) 

The proof of uniform convergence follows from standard spectral theoretic 
arguments [131 Prop. 3.4, p. 12]. □ 
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In the case z = Proposition [S] implies immediately the following theo- 



rem. 



Theorem 3 (Selberg's pre-trace formula). // h satisfies (HI), (H2), 
(H3*), then 

^ 1 

Y^h{pj)\^j{z)\^ ^ — h{p)i&nh{np)pdp+ fc(7^'^)- (179) 

j=o ^ 7e-r-{i} 

which converges absolutely, uniformly in z ^ H^. 

Proof. Use (|114p for the 7 = 1 term. □ 

Using (|169p . the pre-trace formula (|179p becomes 

Y,HPj) = y- / /i(p)tanh(^p)pdp (180) 

+ ^ II II {^3. (5^, gz) - Kgz, gz)] . (181) 

Theorem 4 (Selberg's trace formula). If h satisfies (HI), (H2), (H3*), 
then 

LMft) = ^^y Mp)tanh(.p)pdp+ L 2sinh(r.^/2) ' 

(182) 

which converges absolutely. 

(We will see in the next section (Corollary [1]) that the condition (H3*) 
may in fact be replaced by (H3).) 

Proof. We integrate both sides of the pre-trace formula (|179p over I^\]HI^. By 
the normalization of the eigenfunctions ipj, the left hand side of (|179p 
yields the left hand side of p82p . The first term on the right hand side is 
trivial, and the second term follows from the observation that 

E / fi9-^)dp= I f{z)dfi, (183) 



which allows us to apply Proposition [H] to the inner sum in (|18ip . □ 

Remark 6. The absolute convergence of the sum on the right hand side of 
(|182p only requires (HI), or 

|g(i)| <e~"l*l, Vi>0. (184) 
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One way of seeing this is is that, since g is only evaluated on the discrete 
subset (the length spectrum) 

{l^: 7e^-{l}}cR>o, (185) 

we may replace g by an even C°° (M) function g (for which absolute convergence 
is granted) so that 3(^7) — 5(^7) for all 7, and 

|5(i)| <e~"l*l, Vt>0. (186) 

Remark 7. We may interpret the sum over conjugacy classes in the spirit of 
Propositions [H] and [T] provided h satisfies (HI), (H2), (H3*), we have 

V h{pj) = ^^^^^-^^ I tanh(7rp) p dp + i V [ h{p) n^^ (p) dp. 

(187) 

Alternatively, replace the first term on the right hand side in (I182|) by ()108p . 
then 

^ Area(X) g'jt) ^ ^ gjni,) 

h^''^- ^Jo ^^^W^^,k,tl2sinh(nV2)- ^'''^ 



11 The heat kernel and Weyl's law 

As a first application of the trace formula, we now prove Weyl's law for the 
asymptotic number of eigenvalues Xj below a given A, 

N{X) = #{.? : A, < A}, (189) 

as A 00. 

Proposition 10 (Weyl's law). 

^ Area(A^) 
47r 

Proof. For any /3 > 0, the test function 

h{p) = e'Pp' (191) 
is admissible in the trace formula. The Fourier transform is 

and so ()182p reads in this special case (with Xj ~ Pj + \) 
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^ g.. Arca(M) f°° ^g(^2.U , , , , 
j=o 



71^ 2sinh(n£,/2) ' ^^^^^ 



The sum on the right hand side clearly tends to zero in the limit /3 — > 0. Since 
tanh(7rp) = 1 + 0(e-2'^lpl) for all p e M, we obtain 



OO 



Ee-^^^=^^ + 0(l), P^O. (194) 

The Proposition now follows from a classical Tauberian theorem [18]. □ 

The sum X^jlo ^^'^'^^ represents of course the trace of the heat kernel e^"^. 

Corollary 1. The condition (H3*) in Theorem^ and Remark^ can be re- 
placed by (H3). 

Proof. Weyl's law implies that, for any S > 

OO OO 

+ <oo, i.e., ^(1 + Repj)"2"'^ < OO. (195) 

i=0 ]=Q 

To prove this claim, note that 

OO />oo 

V(l + A,)-i-*/2 ^ / {l + x)-^-^/^dN{x) (196) 

(197) 

+ (^1 + 0^ {l + x)-^-^/^N{x)dx (198) 

(use integration by parts) which is finite since N{x) grows linearly with x. 

If h satisfies (HI), (H2), (H3), then the function h^{p) = h{p)e~'^'' clearly 
satisfies (HI), (H2), (H3*) for any e > 0, with the additional uniform bound 

|/ie(p)|«(l + |Rep|)-^-^ (199) 

where the implied constant is independent of e. By repeating the calculation 
that leads to ([75]) , we obtain the following estimate for the Fourier transform 

of /le, 

|<?,(t)|«e--l*l, (200) 
where the implied constant is again independent of e. Theorem |4] yields 



j=o 

= il + x)-^-^/^N{x) 



x=0 
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(201) 

Due to the above e-uniform bounds, both sides of the trace formula converge 
absolutely, uniformly for all e > 0. We may therefore take the limit e — > 
inside the sums and integral. 



12 The density of closed geodesies 

In the previous section we have used the trace formula to obtain Weyl's law 
on the distribution of eigenvalues A^- . By using the appropriate test function, 
one can similarly work out the asymptotic number of primitive closed geodesic 
with lengths < L, 

n{L) = #{7 eH,:£^< L}. (202) 
In view of Remark [6] we know that, for any S > 0, 

^ £^e-'^^'+''> < 00 (203) 

which implies that, for any e > 0, 

i7(i)<, e^(i+^). (204) 

There is in fact an a priori geometric argument (cf. [13]) which yields this 
bounds with e = 0, but the rough estimate (|204p is sufficient for the following 
argument. 

Let us consider the density of closed geodesies in the interval [a + L,b + L] 
where a and b are fixed and L ^ 00. To avoid technicalities, we will here only 
use smoothed counting functions 



V'l(^) = / i^L{t)dnit) (205) 
■yen, ^0 

where i'Lit) — ipit — L) and ij^ E C^(R). One may think of as a smoothed 
characteristic function of [a,b]. Stronger results for true counting functions 
require a detailed analysis of Selberg's zeta function, which will be introduced 
in Section [141 

Let Pi), ... , pm be those pj with Im pj < 0. The corresponding eigenvalues 
Ao, . . . , Xm are referred to as the small eigenvalues. 

Proposition 11. Lei V S Cg°(M). Then, for L>1, 

tPL{t) dn{t) = I Mt) dn{t) + o , (206) 
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where 

dn(t) = 2J dt. (207) 

Proof. The plan is to apply the trace formula with 

9{t) = ^^^^^ [i^it -L)+ i,{-t - L)] (208) 
which is even and, for L large enough, in C^(M). Hence its Fourier transform, 
h{p) = J J (c(5+'^)* - e(-^+"')*) [i^{t -L) + ^i-t - L)] dt (209) 
= f l(^e^l+'P^*-e^-h+'P)* + e^-l-'P^*-e^h-v)*)^{t-L)dt, (210) 
satisfies (HI), (H2), (H3). Let us begin with the integral 

f i e(^+'^)V(i -L)dt = e^S+'p)^ [ e^^+"'^'^{t) dt. (211) 
Repeated integration by parts yields the upper bound 



I supp i/j 

L{l + \p\~ 



So 

Xr^^"^'^(*-^)'^*«-Z(rW ^'''^ 

This bound is useful for Imp = 0. The other corresponding integrals can be 
estimated in a similar way, to obtain the bounds (assume —1/2 < Imp< 0) 

[ 1 e(-|+'^)V(t - L) dt «M j (215) 

Jr t + \p\) 

f 1 e(-5-'^)V(i - L) dt «;v ,..^^1 |.^ , (216) 
Jr t ^(1 + \P\) 

and 

/J^^^"''''^(*-^)^*<<-L(iTW- ^'''^ 

Therefore, using the above bound with iV = 3, say, yields 
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and 



E ^^Pi) / ^^P) ^^M^P) pdp^OM( — ]. (219) 

The sum of the above terms equals, by the trace formula, the expression 

oo ^ 

E E -^^^7)- (220) 



n 

yeH, n=l 



The a priori bound (|204p tells us that terms with n > 2 (corresponding to 
repetitions of primitive closed geodesies) are of lower order. To be precise, 

oo ^ 

y y -Vl«) «, y - e^L+m+A/n ^^(1+^)72^ (221) 

7eH. n=2 2<n<(L+6)/£„,i„ 

where we assume that ipL is supported in [a + L,b + L], and £min is the length 
of the shortest primitive closed geodesic. Therefore 

y V'l(^) = r Mt) dnit) + O (e^(i+^)/2) . (222) 
The leading order term as L ^ oo is 

E Vl(^) ^ / V'L(i) - dt « — — (223) 
which leads to the improved upper bound for the sum involving repetitions, 

^ 1 ^ e{L+b)/n L/2 

76J?. n=2 2<n<(L+6)/£„i„ 

and hence leads to the desired improved error estimate in (|222p . □ 



13 Trace of the resolvent 

The trace of the resolvent i?(A) = {A + X)^^ is formally 
00 00 

Tri?(A) = E(A - A,)-i = E Mp.), Hp') = ip' - p'Y' , (225) 
j=o j=0 

where p = \ — j as usual. The test function h does not, however, respect 
condition (H3). To overcome this difficulty, we define the regularized resolvent 
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RiX) = (A + \)-^ - {A + K)-^ (226) 
for some fixed A* . The corresponding test function is 

h{p') = {p^-p'Y'-{pl~p'Y\ (227) 
which clearly satisfies (H3), since 



Hp') 



0{p'- 



(228) 



{P^-P'^){pl-P'^) 
We have already encountered the kernel of the regularized resolvent, 

k{z,w) = Gp{z,w)--Gp,{z,w), (229) 

in Section [51 The trace of the regularized resolvent is thus 



Tri?(p)=5][(p^-p/)-i-(pJ-p/ 



(230) 



which, for any fixed p* ^ {±p_,}, is a meromorphic function in C with simple 
poles a,t p — ipj. h is analytic in the strip | Imp'| < a provided a < \ lmp\ < 
I Im p* I where cr > 1/2. 

The trace formula (|187p implies therefore (use formula (|118p for the first 
term on the right hand side, and (|159p for the second) 



2ni 



1=0 



1 



l(/+i) p, -l(Z+i) 



Y E «3-,(p) + C(p,). (231) 



where 



(232) 



converges absolutely, cf. Remark [H] 
Let us rewrite this formula as 



1 1 

^ -r&H, j=o 



9 9 

P - Pj 



9 9 

pi - Pj 



Aiea{M) 
2^ 



E 

(=0 



p-l{l + ^) p, -l(/ + i) 



-C(p,). (233) 



All quantities on the right hand side are meromorphic for all p G C, for every 
fixed p, e C away from the singularities (this is guaranteed for | Imp*| > 1/2). 
Therefore (|233p provides a meromorphic continuation of 



32 Jens Marklof 

nr{p) ^ n^M (234) 

-leH, 

to the whole complex plane. 

Proposition 12. The function nr{p) has a meromorphic continuation to the 
whole complex plane, with 

(i) simple poles at p = po,±pi,±p2, ■ ■ ■ with residue reSp^ nr = and 



(235) 



where pj is the multiplicity of pj . 

(ii) simple poles at p = i[l + ^) with residue 

( Arca(M) + 1) + A. lfl = Q, , ^ 

res,, i,nr = <Ti^ tti ■/ ' 236 

(iii) the functional relation 

nr{p) + nr{-p) ^ ptanh(7rp). (237) 

TT 

Proof, (i) and (ii) are clear, (iii) follows from the identity (|112|) . □ 



14 Selberg's zeta function 

Selberg's zeta function is defined by 

oo 

^(*)= n n (i"^"^''^'^™^ ' ^^^^^ 

which converges absolutely for Re s > 1; this will become clear below, cf. (|24ip 
and ([242)1 . Each factor 

OO 

Jl (l - e-^^(^+")) (239) 

m=0 

converges for all s e C, with zeros at 

s = s^,n = -ni + i{2ii/i^)v, J/SZ, m = 0,1,2,.... (240) 

Note that s^m = ^+iPi/m where Pi^„j are the scattering poles for the hyperbolic 
cylinder 37\EI^- What is more. 



J oo oc _ ^ 



d 

- log 

m— m— 
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with s 



ip, and thus 



Res>l (i.e., Imp < -1/2). (242) 
Recall the the genus is related to the area of M by Area(Al) = 47r((7 — 1). 



Theorem 5. The Selberg zeta Junction can he analytically continued to an 
entire function Z{s) whose zeros are characterized as follows. (We divide the 
set of zeros into two classes, trivial and non-trivialj. 

(i) The non-trivial zeros of Z{s) are located at s — 1 and s — ^ zL ipj {j = 
1, 2,3,.. .) with multiplicity 



2iij if pj = 
ifpj^O. 



(243) 



The zero at s = I (corresponding to j = 0) has multiplicity 1. 
(ii) The trivial zeros are located at at s ~ —I, I = 0, 1, 2, . . . and have multi- 
plicity 2g - 1 for I ^ and 2{g - 1)(2/ + 1) for I > 0. 

Furthermore Z(s) satisfies the functional equation 



Z{s) ^ Z{1 - s) exp 



Proof Equation ((2331) yields 



Mg - 1) 



V tan(7ru) dv 



(244) 



2s- 



1 z^ ' ^ 

j=0 



(s-i)2+p/ 



2(5 -1)E 



9 9 

pi - Pj 



1=0 



S + l / + i + ip* _ 



Note that 



2s- 1 



1 



(i+ipj) S-(i-ipj) 



hence the corresponding residue is 1. Furthermore 

, 2s - 1 
- 2(5 - 1) 



+ C(p,). (245) 



(246) 



s + 1 



(247) 



has residue 2(17— 1)(2Z + 1) at s = —I. Statements (i) and (ii) are now evident. 
The functional relation follows from (|237p . which can be written as 



|-(s) + |-(1 - s) = 47r(g - l)(s - i) tan[^(s - i) 



(248) 
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Integrating this yields 

fS-l/2 



logZ(s) - logZ(l - s) = 47r(g- 1) / w tan(7rw) + c, (249) 

^0 



that is 

s-l/2 



Z(s)/Z(l-s) = exp 



47r(g — 1) / w tan(7ru) + i 

^0 



(250) 



The constant of integration c is determined by setting s = 1/2. Notice that 
the exponential is independent of the path of integration. □ 

One important application of the zeta function is a precise asymptotics 
for the number of primitive closed geodesies of length less than L, L oo; 
we have [T3] 



n{L) = dn{t) + o\^-^j . (251) 

The error estimate is worse than in Proposition [Til since we have replaced the 
smooth test functions by a characteristic function. The asymptotic relation 
(|251|) is often referred to as Prime Geodesic Theorem, due to its similarity 
with the Prime Number Theorem. The proof of ()251|) in fact follows the same 
strategy as in the Prime Number Theorem, where the Selberg zeta function 
plays the role of Riemann's zeta function. 



15 Suggestions for exercises and further reading 

1. Poisson summation. 

a) The Poisson summation formula ([Ij reads in higher dimension d 

f{m) = (252) 

with the Fourier transform 

7(^)= / /(P)e2--Pdp. (253) 



Prove (|252|) for a suitable class of test functions /. 
b) Show that ()252|1 can be written in the form 



Y fim) = Vol(AR') E (254) 



where L is any lattice in and L* its dual lattice. 
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c) Any flat torus can be represented as the quotient L\M'', where the 
Riemannian metric is the usual euclidean metric. Show that the nor- 
malized eigenfunctions of the Laplacian are 

(Pmla;) = Vol(L\M'*)-i/2g2,nm.x (255) 

for every m € L* and work out the corresponding eigenvalues Xj . 

d) Use P54|) to derive a trace formula for 

oo 
3=0 

where pj = y/Xj (this formula is the famous Hardy- Voronoi formula, 

cf. m)- 

Semiclassics. 

a) Show that for p — > oo 



for all fixed r = d{z, w) > 0. Hint: divide the integral ([55]) into the 
ranges [t,2t) and [2t, oo). The second range is easily controlled. For 
the first range, use the Taylor expansion for cosh i at < = r to expand 
the denominator of the integrand. Relation (|256[) can also be obtained 
from the connection of the Legendre function with the confluent hy- 
pergeometric series F{a, 6, c; z), 

^, ^ , ^r{iy+l) e-^-'+i)" ^/l 1 3 1 



r(j.+ 3) (l_e-2r)l/2 \^2'2' 2' 1 



e^ 



(257) 

b) Show that ([256]) is consistent with [Til eq. (41)]. Hint: use the (s,u)- 
coordinates defined in (|138p . 

c) Compare the Gutzwillcr trace formula [11] , [7] , with the Selberg trace 
formula. (Analogues of the ghost of the sphere (Section [6]) for more 
general systems are discussed in |2|.) 

The Riemann-Weil explicit formula. 

a) Compare the Selberg trace formula with the Riemann-Weil explicit 
formula |12[ eq. (6.7)], by identifying Riemann zeros with the square- 



root y Xj — J of eigenvalues Xj of the Laplacian, and logs of prime 

numbers with lengths of closed geodesies. See [3] for more on this, 
b) What is the analogue of the ghost of the sphere? 
Further reading. In this course we have discussed Selberg's trace for- 
mula in the simplest possible set-up, for the spectrum of the Laplacian 
on a compact surface. The full theory, which is only outlined in Selberg's 
original paper [16], is developed in great detail in Hejhal's lecture notes 
|13| . [2], where the following generalizations are discussed. 
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a) The discrete subgroup F may contain elliptic elements, which leads 
to conical singularities on the surface, and reflections (i.e., orientation 
reversing isometrics). Technically more challenging is the treatment of 
groups r which contain parabolic elements. In this case /^\IHI^ is no 
longer compact, and the spectrum has a continuous part, cf. [I4] . 

b) Suppose the Laplacian acts on vector valued functions / : ^ 
which are not invariant under the action of T-y, but satisfy 

T,I = x{i)I yjer (258) 

for some fixed unitary representation x ■ T U(A^). The physical 
interpretation of this set-up, in the case = 1, is that Aharonov- 
Bohm flux lines thread the holes of the surface. 

c) The Laplacian may act on automorphic forms of weight a, which cor- 
responds, in physical terms, to the Hamiltonian for a constant mag- 
netic field B perpendicular to the surface, where the strength of B is 
proportional to a. See Comtet's article [8] for details. 

I also recommend Balazs and Voros' Physics Reports article [1] and the 
books by Buser [4], Iwaniec [15] and Terras [17], which give a beautiful 
introduction to the theory. Readers interested in hyperbolic three-space 
will enjoy the book by Elstrodt, Grunewald and Mennicke 0. Gelfand, 
Graev and Pyatetskii-Shapiro [lOj take a representation-theoretic view on 
Selberg's trace formula. 
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